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Abstract. The New World genus Micrurus contains more than 80 currently recognized species of venomous coral snakes. 
The taxonomy of the South American M. lemniscatus complex is controversial. Within this group, M. lemniscatus, 
M. carvalhoi, M. diutius, M. frontifasciatus, and M. helleri have been treated either as distinct species or subspecies of 
M. lemniscatus. Additional species (M. filiformis, M. isozonus, M. potyguara, M. serranus) have also been proposed to be- 
long to the M. lemniscatus complex but never included in a molecular phylogeny. In the present work, we sequenced four 
mitochondrial (12S, 16S, cyt b, ND4) and one nuclear (Cmos) genes using specimens of M. helleri from the Andean foot- 
hills of Colombia and Peru and M. filiformis from the Colombian Llanos. Supplemented by previously published sequenc- 
es, we inferred the phylogeny of the M. lemniscatus complex using Bayesian and Maximum Likelihood approaches and 
estimated divergence times based on fossil-calibrated nodes. Our results strongly support the monophyly of the M. lemnis- 
catus complex. Furthermore, populations traditionally assigned to M. helleri represent two non-sister lineages, one occur- 
ring along the Andean foothills and the other in lowland Amazonia. As a consequence, we restrict the name M. helleri to 
the populations of the Andean foothills. According to our results, the M. lemniscatus complex diverged from M. surina- 
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mensis during the late Miocene and diversified during the Plio-Pleistocene. 
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Introduction 


The coral snake genus Micrurus WAGLER, 1824 is a spe- 
ciose group of highly venomous snakes belonging to the 
Elapidae family. Currently, more than 80 species are rec- 
ognized that occur mainly in Central and South America, 
with two species encroaching on the southern USA (SILva 
et al. 2016a, WEST et al. 2019, UETZ et al. 2020). Two ma- 
jor groups have been distinguished within Micrurus: The 
first one includes species with a monadal coloration pat- 
tern (color pattern with one black ring between two red 
rings), and the second group comprises mainly species 
with a triadal pattern (three black rings between two red 
rings; ROZE 1996, CAMPBELL & LAMAR 2004, ZAHER et 
al. 2016, JowERs et al. 2019). One additional group, cor- 
responding to the former genus Leptomicrurus SCHMIDT, 


1937, has been included within Micrurus, but the monophy- 
ly of the four species formerly assigned to Leptomicrurus 
has not been examined yet (RENJIFO et al. 2012, ZAHER et 
al. 2016, JowERs et al. 2019). These three groups are sup- 
ported by morphological characters, such as hemipenial 
morphology, relative tail size, and scalation (SLOWINSKI 
1995, ROZE 1996, CAMPBELL & LAMAR 2004). Individual 
Micrurus species are generally diagnosed by scale counts, 
body proportions, and in particular by color pattern (ROZE 
1996, CAMPBELL & LAMAR 2004, SILVA et al. 2016b). How- 
ever, for some groups, these characters have proven to be 
insufficient, both for delimitating species and supraspecific 
groups (ROZE 1996, HARVEY et al. 2003, SILVA et al. 2016b). 

One challenging group within Micrurus is the M. lemnis- 
catus complex. It includes the taxa traditionally identified 
as subspecies of M. lemniscatus, i.e., M. lemniscatus (LIN- 
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NAEUS, 1758), M. carvalhoi ROZE, 1967, M. diutius Bur- 
GER, 1955, M. frontifasciatus (WERNER, 1927), and M. hel- 
leri SCHMIDT & SCHMIDT, 1925 (ROZE 1996, CAMPBELL & 
LAMAR 2004, SILVA et al. 2016b, FLORIANO et al. 2019). This 
group has an entangled taxonomic history. PIRES (2011), 
using morphological data and the largest sampling so far, 
suggested in his unpublished Ph.D. thesis to recognize 
M. carvalhoi, M. diutius, M. frontifasciatus, and M. lemnis- 
catus as full species, and to treat M. helleri as a synonym 
of M. lemniscatus. Additionally, PIRES (2011) suggested the 
existence of an undescribed species, which was later de- 
scribed as M. potyguara PIRES, SILVA, FEITOSA, PRUDENTE, 
PEREIRA-FILHO & ZAHER, 2014. Independently from PIRES 
(2011), the recognition of M. diutius as a full species was 
also proposed by STARACE (2013), based on a morpho- 
logical analysis of populations from French Guiana. Later, 
without explanation, WALLACH et al. (2014) listed M. car- 
valhoi and M. frontifasciatus as full species, and M. diutius 
and M. helleri as junior synonyms of M. lemniscatus. Based 
mainly on molecular evidence, VALENCIA et al. (2016) el- 
evated M. helleri to full species level but kept all other taxa 
as subspecies of M. lemniscatus. Recently, also using mo- 
lecular data, JowERs et al. (2019) concluded that all sub- 
species of M. lemniscatus, perhaps except M. l. frontifascia- 
tus, should be recognized as full species and that there 
is still some cryptic diversity to be described within this 
group. Nevertheless, several recent publications (PIRES et 
al. 2014, SILVA et al. 2016b, TERRIBILE et al. 2018, FLORIANO 
et al. 2019) recognized only M. diutius and M. potyguara, 
besides a polytypic M. lemniscatus with the three subspe- 
cies M. l. carvalhoi, M. I. helleri, and M. l. lemniscatus. As 
a starting point, we will treat all taxa in the M. lemnisca- 
tus complex as full species, recognizing the results of PIRES 
(2011), PIRES et al. (2014), and VALENCIA et al. (2016). 

The Micrurus lemniscatus complex is distributed 
throughout most of the cis-Andean region of South Amer- 
ica, including the eastern Andean foothills, Amazonia, the 
Cerrado and the Atlantic forest, also including the island 
of Trinidad (Fig. 1; ROZE 1996, FEITOSA et al. 2007, PIRES 
et al. 2014, SILVA et al. 2016b, TERRIBILE et al. 2018, Jow- 
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ERS et al. 2019). Initially it was thought that all taxa with- 
in the M. lemniscatus complex had allopatric distributions 
(Figs 1A-B; ROZE 1996, CAMPBELL & LAMAR 2004). How- 
ever, recent data suggest that the ranges of several species 
do overlap. Micrurus diutius and M. lemniscatus seem to 
occur sympatrically in the eastern and western Guiana re- 
gion, as do M. helleri and M. lemniscatus in central Ama- 
zonia, and M. carvalhoi and M. lemniscatus in the Brazilian 
Cerrado (Fig. 1C). In addition, a probable overlap between 
the ranges of M. carvalhoi and M. potyguara was suggested 
for western Brazil (PIRES 2011, PIRES et al. 2014, SILVA et al. 
2016b, TERRIBILE et al. 2018, FLORIANO et al. 2019). 

In addition to the above-mentioned taxa, three further 
species have been proposed to belong to the M. lemnis- 
catus complex by Ferrosa et al. (2007) and TERRIBILE et 
al. (2018): Micrurus filiformis (GUNTHER, 1859), M. iso- 
zonus (COPE, 1860), and M. serranus Harvey, APARI- 
CIO & GONZALEZ, 2003. Micrurus filiformis is distributed 
throughout most of the Amazon Basin and the surround- 
ing Andean foothills, displaying a widely overlapping dis- 
tribution with most other taxa belonging to the M. lemnis- 
catus complex (SCHMIDT & WALKER 1943, ROZE 1996, 
Harvey et al. 2003, CAMPBELL & LAMAR 2004, FEITOSA 
et al. 2007). Micrurus isozonus occurs in the dry forests 
and savannahs of the Guiana Shield, where it could be lo- 
cally sympatric with M. diutius, M. filiformis, M. helleri, 
and M. lemniscatus (ROZE 1996, KOK et al. 2003, CAMP- 
BELL & LAMAR 2004). Without explicit evidence, FEITO- 
sA et al. (2007) suggested that M. isozonus belongs to the 
M. lemniscatus complex, while other authors suggested 
that this species is more closely related to the M. fronta- 
lis complex (ROZE 1996, CAMPBELL & LAMAR 2004). Fi- 
nally, M. serranus occurs in the eastern Andes of Bolivia 
and apparently is not sympatric with any other species of 
the M. lemniscatus complex (HARVEY et al. 2003). TERRI- 
BILE et al. (2018) proposed the inclusion of M. serranus in 
the M. lemniscatus complex, but the original description 
(Harvey et al. 2003) suggested that M. serranus most likely 
belongs to the M. frontalis complex, despite some morpho- 
logical characters resembling the M. lemniscatus complex. 


Figure 1. Proposed distribution ranges of taxa of the Micrurus lemniscatus complex according to (A) Roze (1996), (B) CAMPBELL 
& LAMAR (2004), and (C) FLor1Ano et al. (2019). Blue: M. carvalhoi, red: M. diutius, yellow: M. helleri, green: M. lemniscatus. Red 
circle: putative sympatric occurrence of M. helleri and M. lemniscatus. 
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Relationships within the Micrurus lemniscatus complex 


Until now, molecular phylogenetic analyses included 
only up to four of the six taxa of M. lemniscatus sensu lato 
(M. carvalhoi, M. diutius, M. helleri, M. lemniscatus), and 
none of the additional species suggested to belong to the 
complex (SILVA & SITES 2001, RENJIFO et al. 2012, JOwERS 
et al. 2019). These studies found the studied taxa mono- 
phyletic. However, all analyses were based only on the mi- 
tochondrial ND4 gene, even when more genes were se- 
quenced (e.g., JowERs et al. 2019 sequenced five genes but 
used only ND4). 

In order to better understand the genetic variation with- 
in the M. lemniscatus complex, we present herein a dated 
multilocus analysis. To this end, we used four mitochon- 
drial genes (12S, 16S, ND4, cyt b) and one nuclear locus 
(Cmos). We also included, for the first time, samples of 
M. helleri from the eastern foothills of the Colombian and 
Peruvian Andes and samples of M. filiformis from the Co- 
lombian Llanos. 


Materials and methods 
Sampling 


We studied six samples of Micrurus helleri from the An- 
dean foothills of Colombia, one from the Amazonian 
slopes of Peruvian Andes, and two samples of M. fili- 
formis (Fig. 2; Supplementary Table S1). Tissue samples 
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came from the Banco de ADN y Tejidos de la Biodiversi- 
dad (BTBC) of the Instituto de Genética (IGUN), Univer- 
sidad Nacional de Colombia, corresponding to specimens 
deposited in Colección de Reptiles, Instituto de Ciencias 
Naturales, Universidad Nacional de Colombia (ICN), 
the Instituto Nacional de Salud de Colombia (INSZ), 
Bogota, Colombia, and the Laboratoire Evolution et Di- 
versité Biologique, Université Toulouse, France. We ob- 
tained sequences for five genetic markers, one nuclear lo- 
cus (Cmos) and four mitochondrial genes (125, 16S, cyt b, 
ND4). These genes were chosen to match the dataset of 
JoweRs et al. (2019). Some previously published sequenc- 
es of the triad group were excluded since they appeared 
in exploratory analyses in unexpected positions (i.e., 
deeply divergent from conspecifics) or were unstable be- 
tween analyses, destabilizing branching patterns (Supple- 
mentary Table S2), suggestive of misidentification or sub- 
standard sequence quality. Our final alignment included 
17 specimens belonging to five taxa of the M. lemniscatus 
complex (M. carvalhoi, M. diutius, M. filiformis, M. hel- 
leri, M. lemniscatus) plus 20 additional Micrurus species. 
Since there has been some confusion regarding the iden- 
tity of sequences published by SILVA & SITES (2001), we 
follow RENJIFO et al. (2012) and Jowers et al. (2019), who 
resolved this issue and presented trustworthy taxonomic 
identifications (Supplementary Table S3). For dating pur- 
poses, we used as outgroups sequences of 62 species of 
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Figure 2. Localities of the studied samples of the Micrurus lemniscatus complex and type locality of M. helleri. 
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Table 1. Primers used in this study. * Primers modified by ZAHER et al. (2009) 


Gene Primer Sequence Reference 

128 L1091mod 5 CAAACTAGGATTAGATACCCTACTAT 3’ Kocuer et al. (1989) 
H1557Mop* 5’ GTACRCTTACCWTGTTACGACTT 3’ KNIGHT & MINDELL (1994) 

16S L2510mod (16Sar)* 5’ CCGACTGTTTAMCAAAAACA 3’ PaLuMBI et al. (1991) 
H3056mod (16Sbr)* 5’ CTCCGGTCTGAACTCAGATCACGTRGG 3’ PALuMBI et al. (1991) 
703Botp.mod* 5 TCAAAYATCTCAACCTGATGAAAYTTYGG 3’ Pook et al. (2000) 

gth MVZ16p.mod* 5 GGCAAATAGGAAGTATCAYTCTGGYTT 3’ Pook et al. (2000) 

ND4 ND4ab 5 CACCTATGACTACCAAAAGCTCATGTAGAAGC 3’ — AREVALO et al. (1994) 
H-Leu 5 ATTACTTTTACTTGGATTTGCACCA 3’ STUART & PARHAM (2004) 

bak S77 5 CATGGACTGGGATCAGTTATG 3’ Lawson et al. (2005) 
S78 5 CCTTGGGTGTGATTTTCTCACCT 3’ Lawson et al. (2005) 


the booid and caenophidian radiations, obtained from 
ZAHER et al. (2019). 


Laboratory procedures 


Genomic DNA was extracted using the innuPREP DNA 
Micro Kit (Analytik Jena GmbH, Jena, Germany), follow- 
ing the manufacturer’s protocol. Primers are listed in Table 
1. PCR reactions contained in a final volume of 20 uL 10- 
100 ng of genomic DNA, 1 unit of Taq polymerase (Bioron, 
Ludwigshafen, Germany), 2 uL buffer (as recommended by 
the supplier), 0.5 uM of each primer, and 0.2 mM of each 
dNTP (Fermentas, St. Leon-Rot, Germany). For all genes, 
the same PCR protocol was used. It had an initial dena- 
turation step at 94°C for 5 min, followed by 35 cycles that 
included denaturation at 95°C for 45 s, annealing at 55°C 
for 45 s, and extension at 72°C for 1 min. The cycle ended 
with an extension step at 72°C for 10 min. Purification and 
sequencing protocols followed the conditions described in 
Fritz et al. (2012). 


Molecular analyses 


Sequences were aligned and inspected using GENEIOUS 
9.1.8 (KEARSE et al. 2012) and the implemented MUS- 
CLE algorithm (EDGAR 2004) under default conditions. 
Individual gene alignments were concatenated using SE- 
QUENCE MATRIX 1.8 (Varpya et al. 2011). For phyloge- 
netic analyses, three different partitioning schemes were 
examined using PARTITIONFINDER 2 (LANFEAR et al. 
2017) and the implemented Bayesian Information Criteri- 
on: (1) unpartitioned, (2) partitioned by gene and (3) com- 
plete partition (i.e., partitioned by gene and for protein- 
coding genes by codon position), resulting in the selection 
of the complete partition scheme (Table 2). 

Using the selected partition scheme and the appropri- 
ate evolutionary models, a phylogenetic tree was built with 
the Maximum Likelihood (ML) approach as implemented 
in IQ-TREE 2.0.6 (NGUYEN et al. 2015). Node support was 
estimated using 5000 pseudoreplicates of ultrafast boot- 
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Table 2. Partitions and evolutionary models for the concatenated 
alignment obtained using PartitionFinder 2 (PF) and bModelTest 
(bMT). 


Partition Position PF bMT 
12S 1-540 GTR+4+I+G 121343 
16S 541-1101 GTR+I+G 

cyt b_1 1102-2219\3 GTR+4+I+G 123345 
cyt b2 1103-2219\3 TVM+I+G 

cyt b3 1104-2219\3 TIM+I+G 

ND4_1 2220-3082\3 GTR+4+I+G 

ND4_2 2221-3082\3 TVM+I+G 

ND4_3 2222-3082\3 K81UF+G 

Cmos_1 3083-3652\3 HKY+G 121121 
Cmos_2 3084-3652\3 K80+G 

Cmos_3 3085-3652\3 HKY+G 


strap, considering values of 95% as high support (MINH et 
al. 2013). 

Bayesian Inference (BI) and the Relaxed Molecular 
Clock analysis were run simultaneously in BEAST 2.6.1 
(BOUCKAERT et al. 2019). For BI, we used the best parti- 
tioning scheme from PARTITIONFINDER 2. However, 
evolutionary models were inferred using bMODELTEST 
1.2 (BOUCKAERT & DRUMMOND 2017), exploring all avail- 
able models. In exploratory analyses some parameters were 
stable and effective sample sizes were below 200. Therefore, 
we linked for the software BEAUTIi, included in the BEAST 
package, site and clock models to avoid overparametriza- 
tion. Site models were linked (i) for non-protein-coding 
mitochondrial markers (12S and 16S), (ii) for the protein- 
coding mitochondrial markers (ND4, cyt b), and (iii) the 
protein-coding nuclear marker (Cmos) was treated sepa- 
rately (Table 2). The clock was linked on the one hand for 
the mitochondrial partitions and on the other for the nu- 
clear partition (Cmos). Finally, all partitions were linked 
for the tree model using the Birth-Death Model; the log- 
normal relaxed clock was set as clock model. We used nine 
fossils as node calibration points (Table 3), following Za- 
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Table 3. Fossils used for dating. Dates are given in million years. 


Node Fossil 


Stem Boinae Titanoboa cerrejonensisi 


Stem Colubriformes Procerophis sahnii 
Stem Viperidae Vipera cf. antiqua 
Crotalinae indet. 
Elapid morphotype A 
Incongruelaps iteratus 
Colubridae indet. 


Natricidae incertae sedis 


Stem Crotalinae 
Stem Elapidae 
Stem “Oxyuranine” 
Stem Colubroidea 
Crown Natricidae 


Stem Dipsadidae Paleoheterodon tiheni 


Minimum Maximum Source 


58 93.9 HEAD et al. (2009) 

54 93.9 Race et al. (2008) 

22.1 93.9 SZYNDLAR & BOHME (1993) 
11.2 54 Ivanov (1999) 

24.9 54 McCartney et al. (2014) 
10 54 SCANLON et al. (2003) 

35.2 54 SMITH (2013) 

13.8 54 RAGE & SZYNDLAR (1986) 
12.5 54 Horman (1964, 1977) 


HER et al. (2018, 2019), enforcing monophyly and with log- 
normal distributions with a mean of 1.0 (parameter M) 
and a standard deviation of 1.25 (parameter S). Two inde- 
pendent chains of 100 million generations, sampling eve- 
ry 5000", were run in BEAST, using the CIPRES portal 
(www.phylo.org/; MILLER et al. 2010). Chain convergence 
and burn-in value (10%) were assessed using TRACER 1.7.1 
(RAMBAUT et al. 2018). A consensus tree was summarized 
using TREEANNOTATOR as implemented in BEAST 
2.6.1 (BOUCKAERT et al. 2019). 

Finally, for the individual species uncorrected p distanc- 
es were calculated for the gene with the broadest taxonom- 
ic coverage (ND4) using MEGA X (Kumar et al. 2018) and 
the pairwise deletion option. 


Results 


Both tree building methods revealed the same well-sup- 
ported general topology. Micrurus was found as a highly 
supported clade (Fig. 3), and Sinomicrurus maclellandi was 
identified as sister taxon of Micrurus, albeit with weak sup- 
port (Supplementary Figs Sı and S2). Within Micrurus, 
two deeply divergent and well-supported clades were re- 
covered: One comprised the monad species and the other 
comprised M. narducci plus the triad species. The mono- 
phyly of the triad species clade was only weakly supported. 
The triad clade comprised two well-supported subclades. 
One contained M. dissoleucus, M. mipartitus, M. obscurus, 
and the species of the M. frontalis complex; and the other 
one, M. ortoni, M. surinamensis, and the M. lemniscatus 
complex. 

Both ML and BI approaches revealed the M. lemnisca- 
tus complex as a well-supported clade, and M. surinamen- 
sis as its sister taxon (Figs 3 and 4). All evolutionary rela- 
tionships within the M. lemniscatus complex were identi- 
cal for both analyses and robustly supported (Fig. 3). All 
taxa represented by more than one sequence were recov- 
ered as monophyletic, except for M. helleri. The only cen- 
tral Amazonian sample identified as M. helleri clustered 
with M. carvalhoi, whereas the Colombian and Peruvian 
foothill samples constituted the sister clade of M. diutius, 
rendering M. helleri polyphyletic. 


The M. lemniscatus complex consisted of two well-sup- 
ported clades, one contained M. carvalhoi, M. filiformis, 
M. lemniscatus, and the Amazonian sample of M. helleri; and 
the other comprised M. diutius and the specimens of M. hel- 
leri from the Andean foothills. Within the first clade, M. car- 
valhoi was sister to the Amazonian M. helleri, with M. lem- 
niscatus and M. filiformis as successive sister taxa (Fig. 3). 

Our time-calibrated phylogenetic analysis suggested 
that the M. lemniscatus complex diverged from M. suri- 
namensis in the late Miocene (approx. 5.64 million years 
ago = Mya; Fig. 4; Supplementary Fig. S3). Diversifica- 
tion within the M. lemniscatus complex commenced in the 
Pliocene (approx. 3.86 Mya), while the individual species 
diverged during the Pleistocene, with M. filiformis repre- 
senting the earliest branch-off (approx. 2.33 Mya). The di- 
vergence inferred for the Amazonian M. helleri sample and 
M. carvalhoi was about 1.08 Mya; the split between M. diu- 
tius and the M. helleri samples from the Andean foothills 
was estimated at about 1.78 Mya. 

The lowest uncorrected p distance between Micru- 
rus species for the mitochondrial ND4 gene was between 
M. brasiliensis and M. frontalis (1.58%; Supplementary Ta- 
ble S4). Within the M. lemniscatus complex, the lowest di- 
vergence occurred between the M. helleri samples from the 
Andean foothills and M. diutius (2.41%), and the highest 
between the two clades of M. helleri (8.87%). 


Discussion 
Taxonomy of the Micrurus lemniscatus complex 


Our results suggest that the Micrurus lemniscatus com- 
plex contains at least six species, one of them undescribed. 
Morphology-based research has previously led to a re- 
vised taxonomy of M. lemniscatus sensu lato, but mostly 
focused on whether the individual taxa should be recog- 
nized as species or subspecies (Table 4). In an unpublished 
dissertation, PIRES (2011) proposed to recognize four spe- 
cies: Micrurus carvalhoi, M. diutius, M. frontifasciatus, and 
M. lemniscatus (the latter including M. helleri). However, 
in subsequent publications this proposal was not followed. 
Instead, a more conservative taxonomy was used by sev- 
eral authors, with M. diutius as a full species and a poly- 
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typic M. lemniscatus with the subspecies M. I. carvalhoi, Our results revealed that M. lemniscatus sensu lato is 
M. l. helleri, and M. I. lemniscatus (Pres et al. 2014, Sitva paraphyletic with respect to M. filiformis. Moreover, the 
et al. 2016b, TERRIBILE et al. 2018, FLORIANO et al. 2019). deep divergence between our samples identified as M. hel- 
Other authors (WALLACH et al. 2014, JoweERs et al. 2019) leri from the Andean foothills and from the Amazon Basin 
distinguished a varying number of distinct species within implies that two distinct species are currently lumped to- 
M. lemniscatus sensu lato (Table 4), with M. helleri recog- gether under this name (Figs 3 and 4). One of these species 


nized as another species by JowErs et al. (2019). 
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Figure 3. Maximum Likelihood tree for Micrurus based on 3,618 base pairs of the concatenated alignment of five genes (12S, 16S, 
cyt b, ND4, Cmos). Values above and below branches correspond to ultrafast bootstrap values and Bayesian posterior probabilities. 
Asterisks indicate maximum support under both approaches; dashes that those clades were not recovered in the Bayesian 50% major- 
ity rule consensus tree. Bar colors correspond to Figure 2 and inset photographs. Photos: MaAEL DEwyNTER (M. lemniscatus, French 
Guiana), Lucas Pousa (M. carvalhoi, Brazil); Juan PABLO HurTapo Gomez (M. filiformis, Colombia), ANTOINE Fouquet (M. helleri, 


Peru, AF4455 in the tree; M. diutius, French Guiana). 
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Table 4. Classification of Micrurus lemniscatus sensu lato according to different authors and the present study. 


Pires (2011) 


M. I. carvalhoi 
M. diutius 
M. I. lemniscatus 


M. carvalhoi M. carvalhoi 

M. diutius M. lemniscatus 
M. frontifasciatus M. frontifasciatus 
M. lemniscatus M. I. lemniscatus 


M. I. helleri 


M. lemniscatus 


M. lemniscatus M. lemniscatus 


M. l. lemniscatus 
M. helleri 


WALLACH et al. (2014) SILva et al. (2016b) VALENCIA et al. (2016) JoweRs et al. (2019) This study 


M. l. carvalhoi M. carvalhoi M. carvalhoi 
M. l. diutius M. diutius M. diutius 
M. l. frontifasciatus unclear unclear 


M. lemniscatus M. lemniscatus 


M. helleri M. helleri: two species 


finding that conflicts with the hypothesis of several authors 
(PIRES 2011, SILVA et al. 2016b, TERRIBILE et al. 2018, FLORI- 
ANO et al. 2019) that M. helleri is a subspecies of M. lemnis- 
catus. The second species (from Amazonia) currently iden- 
tified as M. helleri is the sister species of M. carvalhoi. 
Some Micrurus species displaying low genetic distanc- 
es for the mitochondrial ND4 gene (Supplementary Table 
S4), like M. brasiliensis and M. frontalis (1.58%) or M. bal- 
iocoryphus and M. pyrrhocryptus (2.00%), have undergone 
taxonomic revision based on relatively large sampling (S1L- 
VA & SITES 1999). As a result, their species status has been 
widely accepted (HARVEY et al. 2003, SILVA et al. 2016b, 
COSTA & BÉRNILS 2018, NOGUEIRA et al. 2020). The same is 
true for M. fulvius and M. tener (CASTOF et al. 2012, STREI- 
CHER et al. 2016), differing by 2.67%. Some of these diver- 
gences are lower than those among taxa belonging to the 
M. lemniscatus complex (2.41-8.87%; Supplementary Table 
4), supporting the recognition of M. carvalhoi, M. diutius, 
M. filiformis, M. helleri, and M. lemniscatus as full species 


and that an undescribed species exists that has been con- 
fused with M. helleri. 

Micrurus helleri was described from the Andean foot- 
hills of Peru (SCHMIDT & SCHMIDT 1925), with a proposed 
distribution across the Andean foothills of Colombia, Ec- 
uador, Peru, and Bolivia and in westernmost Amazonia 
(Figs 1 and 2; ROZE 1996, CAMPBELL & LAMAR 2004, PIRES 
2011, SILVA et al. 2016b, VALENCIA et al. 2016, TERRIBILE et 
al. 2018, FLORIANO et al. 2019). However, our results sug- 
gest that “M. helleri” is composed of two distinct species. 
Our sample from Peru (AF4455) comes from a locality ap- 
proximately 400 km north of the type locality of M. helleri 
(Fig. 2) but in the same ecoregion (Very Humid Premon- 
tane Forest; BRITTO 2017). Therefore, we restrict the name 
M. helleri to the foothill clade. To our knowledge, there are 
no names available for the Amazonian lineage. Further in- 
vestigations with expanded sampling and additional lines 
of evidence, in particular morphology, should be under- 
taken to properly delimit and name this species. 


: M. surinamensis 
' M. helleri (Foothills) 


: M. diutius 


5.3 


7 M. filiformis 
: M. lemniscatus 
7 “M. helleri” (Amazonia) 


: M. carvalhoi 


0 


Miocene Pliocene 


Pleistocene 


Figure 4. Time tree of the Micrurus lemniscatus complex. Node values are mean dates in millions of years; grey bars indicate 95% 
Highest Posterior Density (HPD) intervals. Inset photo of M. helleri: Juan PABLO HURTADO GÓMEZ. 
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Table 5. Dates for the main clades obtained in the BEAST analysis (95% HPD) and comparison with previous studies. The results 
of JowErs et al. (2019) were based on mutation rates; those of LEE et al. (2016) and ZAHER et al. (2016, 2019) on a fossil-calibrated 
molecular clock. *Approximate dates because exact values were not given by the authors; NR, not recovered. 


Clade Presentstud Jowers etal. LEE et al. ZAHER etal. ZAHER et al. 
y (2019) (2016) (2016) (2019) 

Micrurus 11.79 (9.56-14.09) 31 18:755 19* 14.78 
Monadal clade 5.28 (3.69-7.79) 8.30* 5:3% 4.21 

Triadal clade 8.91 (6.96-10.92) 24 14.51* 14* VT 

M. surinamensis + M. lemniscatus complex 5.64 (4.26-7.29) 14.8* 10.15* NR NR 

M. lemniscatus complex 3.86 (2.82 -5.38) 9.6* NR NR NR 

M. diutius + M. helleri 1.78 (1.04-2.89) NR NR NR NR 


For M. carvalhoi, an additional ND4 sequence is avail- 
able from GenBank (accession number: AF228435, vouch- 
er: IB55598; Supplementary Table S1). It was included in 
exploratory analyses, and M. carvalhoi was found mono- 
phyletic with strong support (Supplementary Fig. S4), in 
agreement with the findings of RENjiFo et al. (2012) but 
contradicting JoweErs et al. (2019), who found M. carvalhoi 
polyphyletic. However, in our trees, the sequence AF228435 
displayed a long branch. Under ML, the support values for 
other clades were lower, and the BI topology changed. We 
assume that the sequence AF228435 shows some base-call- 
ing errors, which is why we did not use it, except for the ex- 
ploratory analysis presented as Supplementary Figure S4. 

Given the mislabeling and misidentification of many 
GenBank sequences associated with the M. lemniscatus 
complex, a few erroneous interpretations about the evolu- 
tion and diversity of the group have been published (e.g., 
cryptic diversity in JowERs et al. 2019; evolution of venom 
composition in SANz et al. 2019). To avoid future misun- 
derstandings, we summarize misidentifications of Gen- 
Bank sequences in Supplementary Table S3. 


Diversification and biogeography in the Micrurus 
lemniscatus complex 


Our fossil calibrated tree (Fig. 4 and Supplementary Fig. 3) 
retrieved ages for the M. lemniscatus complex more recent 
than the ones reported in previous works (Table 5; LEE et 
al. 2016, ZAHER et al. 2016, 2019, JowERs et al. 2019). This 
could be associated with the fact that we analyzed a larger 
number of genes for the ingroup. Nevertheless, our dates 
were closer to the previously published results using fos- 
sil calibrations (LEE et al. 2016, ZAHER et al. 2016, 2019) 
than to the results based on a priori mutation rates (Jow- 
ERS et al. 2019), which advocates for a careful use of priors 
for molecular dating. 

Species in the M. lemniscatus complex are all considered 
semi-aquatic (ROZE 1996, CAMPBELL & LAMAR 2004, FEI- 
TOSA et al. 2007, ALMEIDA et al. 2016, SILVA et al. 2016b). 
These taxa diverged from their fully aquatic sister taxon 
M. surinamensis in the late Miocene (approx. 5.64 Mya; 
Fig. 4). Micrurus surinamensis is co-distributed with the 
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M. lemniscatus complex throughout most of its range, and 
their diversification might have also been driven by eco- 
logical differences. According to our results, the radiation 
of the M. lemniscatus complex began in the Pliocene (ap- 
prox. 3.86 Mya) and could have been triggered by factors 
like changes in river courses (ALBERT et al. 2018, 2020) or 
the Plio-Pleistocene climatic fluctuations and associated 
forest dynamics, as it has been shown for other vertebrates 
(WUsTER et al. 2005, QUIJADA-MASCARENAS et al. 2007, 
Hoorn et al. 2010, VaRGAS-RAMIREZ et al. 2010). 

Further research is needed to clarify the identity and 
distribution of M. helleri sensu stricto and of the unnamed 
Amazonian species. The distribution of M. helleri sensu 
stricto seems to be continuous across the eastern Andean 
foothills of Colombia, Ecuador and Peru (Fig. 2), as it has 
been found for other snakes (CAMPBELL & LAMAR 2004, 
ARTEAGA et al. 2018), for instance Bothrocophias microph- 
thalmus (Copre, 1875), Dipsas peruana (BOETTGER, 1898) 
sensu lato or D. vermiculata (BOULENGER, 1885), and sev- 
eral other animal groups endemic to this region (LYNCH et 
al. 1997, PATTERSON et al. 1998, KaTTAN et al. 2004, RON et 
al. 2011). Additionally, there are records further north in 
Venezuela, in the Merida foothills, identified as M. lemnis- 
catus (BARRIO-AMOROS & CALCANO 2003, NATERA-Mu- 
Maw et al. 2015). These records could also refer to M. hel- 
leri sensu stricto. The distribution of the unnamed Amazo- 
nian species is unclear but could be identical to the whole 
range originally assigned to M. helleri in central and west- 
ern Amazonia (Fig. 1). 
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based on 3,618 base pairs of the concatenated alignment of five 
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Supplementary Figure S2. Complete Bayesian tree based on 3,618 
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cyt b, ND4, Cmos). 


Supplementary Figure S3. Complete Bayesian time tree with 
mean dates above nodes in million years before present. 


Supplementary Figure S4. Complete Maximum Likelihood tree 
including the specimen Micrurus carvalhoi 1B55598. 
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